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We study the ground state of hard-core bosons with nearest-neighbor hopping and nearest- 
' neighbor interactions on the triangular and Kagome lattices by mapping to a system of spins (S = i), 

, which we analyze using spin-wave theory. We find that the both lattices display superfluid and 

■ supersolid (a coexistence of superfluid and solid) order as the parameters and filling are varied. 
t-H ' Quantum fluctuations seem large enough in the Kagome system to raise the interesting possibility 

I . of a disordered ground state. 

3 ■ 

hj ; 

O ■ INTRODUCTION 

m ; 

A supersolid is a state of matter which simultaneously exhibits both solid and superfluid properties. That is to 
, say, it displays both long-ranged positional order as well as finite superfluid density and, naively, off-diagonal long- 

■ ranged order (ODLRO). The intriguing suggestion by Andreev and Lifshitz [jl| that vacancies in solid 4 He might 
Bose condense in the vicinity of the melting line has, to our knowledge, never been experimentally verified 
Nonetheless, a sizeable literature has developed on the theoretical properties of supersolids ^HUl- I n two dimensions, 
the physics of Josephson Junction Arrays |l7]] has also stimulated the theoretical study of supersolids. Most of the 
work is based on the contributions of Matsuda and Tsuneto |l8| and of Liu and Fisher |l9| , who established some key 

Q\ \ concepts in the theory of lattice-based supersolid models. Of central importance is the mapping between a hard-core 
lattice Bose gas and a spin-i quantum magnet: 



o 
o 



4 <-> St a t Si n t = a\ ai <-> S- - \ . (1) 



e. 

Thus, an occupied site is represented by an up spin, while an empty site is represented by a down spin. An interacting 
hard core lattice Bose gas with nearest neighbor hopping t, nearest neighbor repulsion V, and chemical potential fi is 
thereby equivalent to the anisotropic S — h Heisenberg model 



where Jy = V is an antiferromagnetic longitudinal exchange, J± = —It is a ferromagnetic transverse exchange, and 
is an external magnetic field (z is the lattice coordination number). The spin model exhibits a global 
U(l) symmetry with respect to rotations about the z axis, which of course means total particle number conservation 
in the boson language. The boson condensate order parameter is related to the transverse magnetization density via 
(a|) = (S^), while the boson compressibility K = dn/dfj, is the magnetic susceptibility x = dM z /dH. Liu and Fisher 
identified four phases of interest: (i) a normal fluid, in which the magnetization is uniform and in the z direction, 
(ii) a normal solid, in which the magnetization lies along z yet is spatially modulated at some wavevector k, (iii) 
a superfluid, in which the magnetization is uniform and has a component which lies in the x-y plane, and (iv) a 
supersolid, in which there simultaneously exists a nonzero transverse component to the magnetization M±, as well 
as a spatial modulation of the longitudinal magnetization M z . (An incompressible normal fluid is also called a Mott 
insulator.) 

If one relaxes the hard core constraint in favor of a finite on-site repulsion U , one obtains the Bose Hubbard model 

H = -ty^ j (a\a i + ajaj - fj,y^/n,j + ~U ^^{nj - 1) + V njUj . (3) 
(ij) » • <y> 

This model has been extensively studied since the seminal work of Fisher, Weichman, Grinstein, and Fisher poj , who 
considered the model with V = in the context of a superconductor-insulator transition. A study of this model in 
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the presence of disorder has led to an understanding of the Bose glass pi]] . On a two-dimensional square lattice, 
and for V ^ 0, the model was studied at T = 0, for both finite and infinite U, by Scalettar et al O], Bruder e£ 
fl, and van Otterlo et al using mean field theory and quantum Monte Carlo techniques. To summarize their 
results, no supersolid phase is observed at half filling ((n) = where a first order transition occurs as a function 
of V between superfluid (small V) and a Neel solid (large V, k = (tt,tt)). (Large next-nearest neighbor repulsion V 
stabilizes a striped phase, the colinear solid.) Away from half filling, there is no normal solid phase, and the transition 
is instead from superfluid to supersolid. The supersolid phase exhibits both a peak in the static structure factor S(k) 
at the Neel vector (and properly proportional to the lattice volume), and a nonzero value of the superfluid density 
p s . Again, next nearest neighbor V' can stabilize a striped supersolid phase with anisotropic p s . One can also obtain 
Mott insulating phases with fractional filling in the presence of next nearest neighbor interactions. 

In this paper, we will investigate the properties of the model in eq. |] on frustrated two-dimensional lattices. We 
are motivated by the fascinating interplay between frustration, quantum fluctuations, order, and disorder which has 
been seen in quantum magnetism. 

Frustration enhances the effects of quantum fluctuations. Indeed, as early as 1973, Fazekas and Anderson p3 24 
raised the possibility that for such systems, quantum fluctuations might destroy long-ranged antiferromagnetic order 
even at zero temperature. In many cases, frustration leads to an infinite degeneracy at the classical (or mean field) 
level not associated with any continuous symmetry of the Hamiltonian itself. In these cases, it is left to quantum 
(or thermal) fluctuations to lift this degeneracy and select a unique ground state | |25| , |26| , sometimes with long-ranged 
order. Our models exhibit both a depletion (but not unambiguous destruction) of order due to quantum fluctuations, 
as well as the phenomenon of "order by disorder" . 

In our work, we will choose the units of energy to be Jm, writing A = t/V = J±/2Ju, and h = H/J±. We will be 
following closely the analysis of the anisotropic triangular lattice antiferromagnet by Kleine, Muller-Hartmann, Frahm, 
and Fazekas (KMFF) | p7| , who performed a mean field (S = oo limit) and spin wave theory (order 1/5 corrections to 
mean field) analysis. Comtemporaneously with KMFF, Chubukov and Golosov [ p8[ derived the spin-wave expansion 
for an isotropic Heisenberg antiferromagnet in a magnetic field, while Sheng and Henley ]^9| obtained the spin-wave 
theory for the anisotropic antiferromagnet in the absence of a field. 

The mean field phase diagram is shown in figure 1 (both the triangular and Kagome lattices have the same mean 
field phase diagram up to a rescaling of h). Notice that the supersolid phase appears in a broad region of A and 
filling. The reason the supersolid is so robust is that the lattice frustrates a full condensation into a solid. Generically, 
frustrated lattices might be good places to look for this phase. 

Let us briefly concentrate on h = before describing the entire phase diagram. We will be assuming a three 
sublattice structure throughout. The mean field state is then described by three polar and three azimuthal angles: 
(#a; @h,9c, <j>Ai 0b> <fic), an d is invariant under uniform rotation of the azimuths. 

Due to the ferromagnetic coupling in the x — y spin directions the mean field solution is always coplanar. Just 
as in KMFF, there is a one-parameter family of degenerate mean field solutions in the zero-field case (originally 
found by Miyashita and Kawamura pc| ). The A sublattice polar angle 9 A may be chosen as the free parameter; 
spin wave theory (SWT) is necessary to lift the degeneracy and uncover the true ground state. Figure 2 shows the 
ground state energy in SWT as a function of 8 A for the triangular lattice at A = 0.25. Using SWT we also compute 
the fluctuations of the spins, and the consequent quantum-corrected magnetization and the solid and ODLRO order 
parameters. Figure 3 illustrates these quantities in mean field and to leading order in SWT (where S has been set 
equal to |) as a function of A for the triangular lattice. It is clear that the quantum corrected S z is very close to 
zero for all A, reflecting the fact that at h — the lattice is half-filled. Two sublattices acquire large corrections due 
to quantum fluctuations (even in the Ising limit A — > 0), while the third has only small quantum corrections. This is 
very similar to the fully antiferromagnetic case studied by KMFF. Therefore, even at S = h the solid order survives. 
The off-diagonal order parameter S x is reduced in magnitude by quantum corrections, but goes to zero only as A— >0. 
Supersolid order survives quantum fluctuations for the triangular lattice, at least in this order of SWT. 

Let us now consider the Kagome lattice. There is a qualitative difference between the antiferro-ferromagnetic case 
considered here, A > 0, and the fully antiferromagnetic case A < 0, which has been exhaustively explored for the 
Heisenberg limit A = — 1 (for a partial set of references, see [^l]-^8| ) . For the fully antiferromagnetic case there are 
local motions of the spins that move the system on the degneracy submanifold, leading to a much larger ground state 
degeneracy for the Kagome lattice than for the triangular lattice. However, for A > 0, the ferromagnetic transverse 
interaction eliminates the possibility of these local motions, resulting in a ground state degeneracy parametrized only 
by 9 A , just as in the triangular lattice. 

We carried out SWT for the two long-range ordered configurations shown in figure 4 - a three sublattice "<f = 0" 
state, and a nine sublattice \/3 x \/3 structure B3,B9], respectively. These states have the same mean field ground 



2 



state energy. It will turn out that two of the three sublattices have the same spin orientation in the ground state 
for both lattices. With this proviso, note that the q = structure has one two-fold axis and a mirror plane (point 
group C-2 V ), while the a/3 x a/3 structure has a six- fold axis and a mirror plane (point group Cq v ). Once again SWT 
selects the true ground state. When quantum fluctuations are accounted for, we find that the a/3 x a/3 structure 
always has lower energy than the q — structure, to the numerical accuracy of our calculations. More importantly, 
the fluctuations of the spins on six of the nine sublattices diverge in the limit h^O, as shown in figure 5. This 
divergence is the consequence of a flat (dispersionless) mode at zero energy as h— >0. Higher order terms in the spin 
wave expansion will lift this mode and remove the divergence 25 3J|. However, the fluctuations remain large for 
S = \ , as we estimate in section III. This indicates that quantum flucuations may be strong enough to wash out any 
order, including ODLRO, on two of the sublattices, which raises the intriguing possibility of a partially disordered 
ground state at T = Q. It must be emphasized that we have not demonstrated that this is so: the fluctuations could 
be correlated between different sites, there could be long-range order with a larger unit cell, a condensed array of 
vortices, etc. This problem merits further study, with e.g. quantum Monte Carlo methods. 

Let us now turn to a fuller description of figure 1, where four types of mean field states are present (solid lines 
indicate first order transitions). We adopt the nomenclature of ref. (see fig. 12 of this reference for comparison). 
At high fields h the system is in a Mott phase - incompressible and fully polarized. As the field is lowered, for any 
A > 0, the system enters a compressible superfluid phase (SF), with 9 A = 9 B = 9 C > 0. A first order transition from 
the superfluid to an incompressible Neel solid (NS) at filling fraction | (magnetization per site M z — ±\S) occurs 
for A < ~. Finally, the supersolid phase (SS) exists for A < i between the two symmetry- related Neel solid lobes. 
There is a tricritical point at (A*, h*) — (^,3). Quantum fluctuations will modify these mean field phase boundaries. 
Since the superfluid state, we have found, benefits the most from spin wave energy corrections, it will encroach on its 
neighbors as S decreases from oo. Increasing h tends to suppress quantum fluctuations. 

We will present each of the above results in more detail in the rest of this paper. Section I concentrates on the 
mean field theory and the mean field phase diagram. Section II describes the selection of the true ground state by 
quantum fluctuations at h = 0, and the form of the spin-wave excitations for arbitrary h. Section III presents the 
suppression of order by quantum fluctuations. We end with our conclusions, connections to experimental work, and 
open questions in Section IV. 



I. MEAN FIELD THEORY 



The mean field limit is obtained by setting S = oo. Formally we first generalize the model from S = ^ to a model 
with a spin-S 1 at each site. We then represent each spin as a classical vector of magnitude S, S a — Sft a , where Cl is 
a unit vector in three dimensions. We rescale the magnetic field by S and write 

n MF /s 2 = Yl n i QZ i - A E(^ fi i + n ? n ?) - h J2 n i ■ w 

(ij) (ij) i 

In the mean field solution, all spins lie in the x-z plane. Furthermore, it turns out that the triangular and both 
Kagome structures have the same mean field energy, to within a constant factor (up to a separate rescaling of h in 
the case of the Kagome structures) , so we will consider all three cases simultaneously. Note that coplanar states have 
been selected at the mean field level for the Kagome lattice, in contrast to the fully antiferromagnetic case. 



A. Zero Field 



Consider first mean field ground states on the triangular and Kagome lattices that have a three-sublattice structure. 
The three specific cases are the usual sublattice structure on the triangular lattice, the q = structure on the Kagome 
lattice, and the a/3x a/3 structure on the Kagome lattice. The nine sublattices of the a/3 x a/3 structure are organized 
into three groups (A,B,C) of three, so that an A site has two B and two C neighbors. Thus, the energy per site, in 
units of S 2 , is 

e MF = E MF /NS 2 = (cos 6 A cos 9 B + cos 9 B cos 8 C + cos 9 C cos 9 A ) 



3 



— A(sin6* A sin# B + sin6> B sin6> c + sin 0c sin # A ) (5) 



on the triangular lattice and z Kas / z Tri = | this value on the Kagome lattice. 



Miyashita and Kawamura 1 30 1 have shown that there is a one-parameter family of degenerate ground states for this 
classical problem for arbitrary A, which does not seem to be related to any obvious symmetry of the model. Following 
KMFF, and writing = 0i — (9 A + 9 B + 9 C ), and defining the two-dimensional vectors /2j = (sin/3i, cos/3i), we can 
write the mean field energy in terms of the two-component vector jl = fi A + fI B + /2 C : 

e MF = 1(1- A)(/! 2 -3) + i(l + A)/i y . (6) 

Therefore e M F, while nominally depending on the three angles 9i, actually depends only on two combinations of them, 
leaving one parameter free. We can then parameterize the degenerate ground states by 9 Al by defining 9 B — 6 — 
and 9 C = e + S, where 

tan # 4 

tane = 



A 

—A cos 9 » 



^-(l-A^cos 2 9 A 

cos 6 = = (7) 

(1-A),/1-(1- A 2 ) cos 2 9 a 

It is easy to verify that 9 A can lie in the range cos -1 ( i~5^" i-A^ ) — ^ A — 

As A increases from zero, the range of possible 9 A is compressed, and the differences \9 B — 9 A \, \9 C — 9 A \ shrink, until 
at A = i the ground state is colinear with 9 A = 9 B = 9 C = ^ir - a featureless superfluid. For A < i, in the case of 
zero field, the true ground state will be selected by quantum fluctuations. Minimizing e MF gives fj, x = (A+ 1)/(A — 1), 
jj, y = 0, and 



- A + A^ 
1 - A 



(8) 



B. Non-zero Field 



Turning on a field h adds an energy 

Ae MP = — |/i(cos# A + cos0 B + cos# c ) (9) 

per site and lifts the degeneracy described in the previous subsection, producing a unique mean field ground state 
(unlike in the Heisenberg antiferromagnet |Q ) . We find that minimization generally leads to a state where (without 
loss of generality) 9 A — 9 B ^ 9 C in the supersolid phase. The results are plotted in figure 6. A rescaling of the field for 
the Kagome lattice (h Kas = |/i T ri) makes the entire mean field phase diagram identical, and we have therefore shown 
only the triangular lattice results. 

The phase diagram has already been shown in figure 1. Let us keep at a particular value of A and turn up the 
field h. At zero field there are two regimes, the superfluid with no solid order for A > |, and the supersolid with 
both solid and ODLRO (T = 0) for A < |. For A > |, the ground state remains a uniform superfluid, though M z 
becomes nonzero, as h is increased. The spins cant at an angle 9 = cos _1 (/i/6(l + A)). The energy in this phase is 

Eventually, for h > 6(1 + A), every site has the maximum possible S z , and the system is in the MI phase, with 
each site fully occupied with one boson. Borrowing from spin-wave results derived in section II, the linear instability 
of the SF phase occurs at 



hn =6x/(l + A)(l-2A) (11) 
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Next focus on a specific A < |, and increase the field h from zero. At zero field, the one-parameter degeneracy of 
the mean field ground states has to be lifted by quantum fluctuations, which select a particular 9 A . However, it turns 
out that we can recover this 9 A by considering nonzero h and taking the limit h— >0, which gives 

f — 2A 

cos 2 # A = T -^. (12) 

It seems surprising that the ground state selected by quantum fluctuations can be predicted by an entirely classical 
calculation. A plausible (though non-rigorous) argument will be provided for this in terms of spin-wave theory in the 
next section. 

As h is increased from zero, 9 A and 9 C change (recall 8 A — 9 B throughout). At a certain critical field ft. c i(A), Sc 
points exactly along the — z direction (9 C = 7r), while Sa,b point along the z direction (0 A = 9 B = 0). This critical field 
h c i can be analytically determined by the following consideration: the point 9 A = 0, 9 C — ir is always a stationary 
point of the mean field energy. However, for h < h c \ it is a saddle point, while for h > h c \ it is a minimum. Therefore 
the second derivative matrix of e MF (# A , 9 C ) should have a zero eigenvalue at h = h c x- Setting the determinant to zero 
gives 



h c i(A) = |(2 + A - - 4A - 7A 2 ) (13) 

For h < h c i(A), we are in the supersolid phase. Just above h c i, however, the system is exactly | filled, and its 
magnetic susceptibility is zero (the compressibility of the corrresponding boson system is zero). There is no superfluid 
order and the system is again a Mott insulator, the Neel solid . It will be seen in section III that the exact filling of 
| survives quantum fluctuations. Since the superfluid order goes continuously to zero below the transition, it is clear 
that this transition is second-order within in mean field theory. The energy of the Neel solid phase is 

e£ F sl =-l-^ (14) 



One can furthermore determine that the Neel solid phase is linearly stable for |(2 + A — \/4 — 4A — 7A 2 ) < h < 
|(2 + A + V4-4A-7A 2 ). 
Further increasing h, we find a critical field 



h c2 (A) = 2(1 + A) + 4^/(1 + A)(l - 2A) (15) 

beyond which the canted superfluid becomes energetically favored. The transition is first-order since it is far from any 
linear instabilities. Finally, a second order line at h > ft, C 3(A) = 6(1 + A), signals the boundary between superfluid 
and fully polarized Mott phases. 

When A = i, we can solve analytically to find 

cos# A = \h 
cos8 c = —\h 

This is exactly the same as the energy of the SF phase at A = | (which is characterized by all the angles satisfying 
co$>9 = \K), which marks this vertical line as a first order line. Note the tricritical point at A = |, h — 3, where two 
first order lines (with infinite slope) meet a second-order line (with finite slope). 

The situation is completely symmetric with respect to the sign of h, with the Neel solid phase now existing at 
i filling for h < 0. At the mean field level, the A-axis is a first-order transition line up to A = i, since S z is 
discontinuous across it. However, for S = 4 it apparently becomes continuous, at least for the triangular lattice. 

The phase diagram has some similarities to the classic picture of the Mott lobes surrounded by superfluid described 
by Fisher, Weichman, Grinstein, and Fisher (FWGF) p0| . However, there are important differences. FWGF con- 
sidered local Hubbard repulsion, whereas our extended Bose Hubbard model we consider affords the possibility of 
incompressible Mott phases at fillings 0, |, |, and 1. The Neel solid phase found by Scalettar et al , for example, exists 
at filling i . Fractional fillings have also been observed in the square lattice with frustrating longer range interactions 
in ref. ||. Also, the transitions ^from the fractional filling MI phases to the (canted) superfluid are first order. Finally, 
and most notably, the entire region between the two fractional Mott lobes is taken over by the supersolid, and the 
supersolid gives way to the superfluid only beyond a hopping t > \V ■ 
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II. SPIN- WAVE THEORY 



We now develop the spin wave theory (SWT) for this problem. If h / 0, there is a unique ground state (up to 
permutations of the sublattices) . When h = 0, the ground state manifold is parameterized by 9 A , and the other angles, 
9 B and 9 C (in general not equal), can be determined from 9 A and A using eq. 0. We implement SWT in the usual 
way: by first performing local rotations of the spins so that the mean field directions point along the local z-axis. 
Since all the spins are assumed to lie on the x-z plane, we can do this by a rotation about the y-axis. Labelling the 
local frame spins with a tilde, we have 

S% =caB0„S% + sin V S Z S 

Rv R.v R.v 

si = si 

Rv Rv 

S% =-sm6v§% + cos V S% (17) 

Rv Rv Rv v ' 

where the subscript R labels a Bravais lattice site, v a basis element, and 9 V is 9 a.b.c depending on the mean field 
orientation of the v sublattice. The triangular lattice and the q = structure on the Kagome lattice have three 
sublattices while the a/3 x a/3 structure on the Kagome lattice has nine. 
We now describe the spin operators in terms of Holstein-Primakoff bosons 



St = i/A \/2S-ibt %b- 

Rv Y Rv V r Rv Y Rv 

= V2SV4 +o(s- 1/2 ) 
t „/, „/,t 

Ri 

2S ^ nv + 0(S- 1/2 ) 
§1 = ihl ip . - S . (18) 

Rv ^Rv^Rv v ' 



SZ = 4 2S~tbl ib~ tbl 

Rv V y Rv y Rv Y Rv 



The Hamiltonian (restoring h for generality) is now written in Fourier space as 



n svf = E + ^Sj2^Hk)(t *)H>(k):+0(S°) (19) 



k 

(note the normal ordering) where the energy Eq is given by 



E = (N/K)S 2 
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(20) 



Here, N is the total number of lattice sites, K is the number of sublattices, and we define the quantities 



X uv < = cos 9 V cos 9 V > — A sin 9 V sin 9 V < 
Y vv i = sin 9 V sin 6 U > — A cos 9 V cos 9 U , 
h v = hcos9 u (21) 



and z vv i is the number of iS sublattice neighbors each v sublattice site has. The vector k lives in the first Brillouin 
zone of the reciprocal lattice, and 

¥(k) = (4(k),4(k),...,i J Uk),M-h^2(-k),..-,^K(-k) ) ■ (22) 

The matrices M and N have diagonal and off-diagonal elements given by 

v' 

M vv ,{k) = \{Y VV , - A)f uv ,(k) 
N vv {k) = 

Nw>{k) = \{Yw' + &)fw'(k) (23) 
where the function /„„/ (k) is given by the following sum, 
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U'(k) =5^'exp(**.^ (24) 

5 

where the prime on the sum indicates that the sum is over nearest neighbor vectors connecting a v sublattice site to 
a v' sublattice site. 

Now we perform a Bogoliubov transformation, which amounts to finding a rank 2K matrix T satisfying AT = A, 
with 

\VKxK —i-KxK ) 

as well as AT A7i sw T = u), a nonnegative diagonal matrix with identical upper left and lower right blocks. The 
ijj v {k) are the spin wave frequencies. The spin- wave correction to the ground state energy is given by 

K 

AE = ^EE( W "W _ M vv {k)) . (26) 



A. Ground State Selection 



Figure 2 shows the ground state energy, including spin wave corrections, for the triangular lattice at h = as a 
function of 9 A . Since the classical energy is independent of 9 A all the variation comes from the SW correction. It is 
clear that there are two possible values of 9 A which minimize the ground state energy, one of them lying at the edge of 
the allowed range of 9 A . However, the two minima turn out to be physically identical, and correspond to a relabeling 
of the sublattices. We call the minimizing value of 9 A which lies away from the edge of its allowed range 0* . 

The value of 9 A can be determined by purely classical arguments, by extremizing the value of S z oc 9 A + (9 B + 9 C 
along the degeneracy submanifold. We have checked that this is so by a comparison of the analytic expression 
9* A = cos -1 — 2A)/(1 — A 2 ) with the ground state energy curves obtained from SWT. We now provide an 
argument indicating why this might be the case. 

Our argument will be the following: For a generic 9 A there is only one mode of zero energy at k = 0, whereas if 
= there are two zero modes (to spin- wave order). The ground state energy is the sum of the energies of all the 
modes, and does not depend on just the zero modes. However, we think it plausible that having more modes of zero 
energy at k — drags down the energies of all the k modes, thus reducing the full ground state energy. Of course this 
argument is not rigorous |4C| ] , and there may be counterexamples that we are not aware of. 

Let us go on to show the first statement about the number of modes of zero energy. We only sketch the argument 
here: We treat the issue with more generality and greater detail in the appendix. It is helpful to think in terms of 
the coherent states path integral p2] | 

D[6 A ,6 s ,6 a ,<i> A ,<f> B ,<t>c\ e~ A 

(27) 



A = J dr [y^iS cos O^t^^t) + H[0,> 



where TC stands for the spin hamiltonian of equation [| written in terms of 9 and <fi. The first term is the Berry phase 
contribution to the path integral, which also makes Scos9 the momentum canonically conjugate to <j). 

Let us concentrate on just the k = modes. There are three (f> variables and three conjugate 9 variables. We choose 
to call one of the <p variables 4>o corresponding to an overall rotation of all the spins around the z-axis, and call the 
remaining (j)s <f>\ and <p2- 

Choose any particular ground state labelled by 9 A . The hamiltonian for small deviations from the ground state 
configuration is now given by 

U = \9 T M B 9 + \<fM^ (28) 

where 9 T = (59 A , 59 B , S9 C ) and 4> T = (<f>o, <fii, 02 )• The first row and column of are zero since 4>q does not appear 
in the hamiltonian. 

The general procedure for finding the normal modes is the following: 
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• Diagonalize the 2x2 block of M<p and rescale the resulting eigenvectors so that the rescaled M'^ becomes a unit 
matrix in the 2x2 block. Call the rescaled <j> variables ip±. 

• Use the Berry's phase terms to identify the canonically conjugate momenta to <fro, 4>± as Po, P± respectively. 

• Reexpress the matrix Mg as a matrix Mp. 

• Since 0o is cyclic, its canonically conjugate momentum Pq is conserved. Treat Pq as a constant and form linear 
combinations P± = P± + a±P$ such that the off-diagonal terms containing Pq are eliminated. A diagonal term 
multiplying Pq remains. 

• Now diagonalize the lower 2x2 block of Mp. The two eigenvalues of Mp are the squares of the energies of the 
normal modes of the hamiltonian. The third normal mode corresponds to a uniform rotation of all the spins 
around the z-axis coupled to a change in S z , and its energy is always zero regardless of the coefficient of Pq . 

The above is true even if there is no ground state degeneracy. If there is ground state degeneracy in the 9 subspace 
the matrix Mp has a null eigenvector vq. For generic 9 A , Pq has nonzero overlap with vq. In this case the above 
procedure always produces a zero coefficient for P§. Thus there is still only one mode of zero energy, the S z mode. 
Another way of seeing this is to recognize that as long as Pq and Vq have nonzero overlap, one can always rescale vq 
so that it becomes canonically conjugate to <j> . One then has to modify to keep them independent of vq in the 
Poisson Bracket sense. 

However, if P and vo have zero overlap, the null vector must lie in the subspace of P± . This means one of the 
eigenvalues of the lower 2x2 block of Mp must be zero, implying that there is another zero mode of the hamiltonian, 
apart from the S z mode. The condition for Pq and vq to have zero overlap is identical to = along the degeneracy 
direction, which is the same as extremizing S z . 

In brief, if ||j ^ 0, there is only one zero mode, whereas if ||j = 0, there are two zero modes at k = 0. A rigorous 
derivation is supplied in the Appendix. 

We have done explicit calculations to verify all these statements for the triangular lattice. If one computes the 
frequencies of the k = modes at a non-optimal A , one finds one mode of zero energy (the S z mode), and two modes 
of non-zero energy, neither of which is exactly along the degeneracy direction. However, at precisely the optimal 
points, there appear two modes with zero energy, one of which is the S z mode and the other exactly the degeneracy 
mode. The third mode still has nonzero energy. 

In particular, the argument makes no assumptions about the interactions other than that they should conserve 
total S z . So this result should hold even for site dilution or longer range interactions. Of course, for site dilution, one 
should focus on a particular realization of randomness and look at S z over the degeneracy subspace. In general, our 
argument is that the search for the true ground state can be restricted to the points on the degeneracy submanifold where 
all conserved quantities commute in the Poisson Bracket sense with the generators of motion along the degeneracy 
submanifold. 

This criterion can fail if the k ^ modes do not follow the behavior of the k = modes. Also, if there are more 
conserved quantities than degeneracy directions, not all of them may be extremized at the true ground state. 
With 9 A and 6 B taking on the value quoted above, we find that 9 C is at the extreme edge of its allowed range: 

cos^-^/^f (29) 



B. Spin- Wave Dispersions 



1. Triangular lattice 

The BZ of the triangular lattice is shown in figure 7, with the lattice and reciprocal lattice vectors being 

ei =a(l,0) 
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G 2 = i^(0,l) (30) 
avo 

where a is the lattice spacing. 

As usual, we first concentrate on h = 0. Figure 8 shows the spin- wave dispersion for A > i. Since there is only 
one sublattice, there is only one mode. However, plotting it on the reduced BZ of the sublattice problem forces us to 
fold it back and represent it as three modes. In this scheme, there is one gapless mode, which is the Goldstone mode 
coresponding to the density fluctuations of the bosons. However, as A approaches |, a "roton" minimum develops, 
precisely at the wave- vectors corresponding to the sublattice structure, as demonstrated in figure 9. At exactly A = i, 
this becomes a gapless mode, heralding the transition to the supersolid phase. 

Now consider the dispersion at h = in the SS phase, which is shown in figure 10. There are two gapless modes 
within SWT. One is the standard density fluctuation, whereas the other corresponds to the degeneracy mode. The 
degeneracy mode will be shown to acquire a gap to higher order in 1/5 in the next subsection p5|j39| ]. There is a 
third "optical" mode to complete the count of the sublattice degrees of freedom. The energy scale of the two low-lying 
modes is A, while the optical mode has an energy scale of 1. As we move towards A — the two low-lying modes get 
softer, until they become completely flat at A = 0. 

However, when we turn on a field the degeneracy mode becomes gapped, as shown in figure 11. Let us now 
investigate the modes at nonzero field in the other phases, in particular at the transitions. The SWT for the canted 
spin phase can be analysed analytically, since there is only one sublattice. We find that the SW dispersion is 



u(k) = zsJ A(l - 7fc -)(A + (1 - * W-) 



36(1 + A)' 



Z' 
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where z = 6 for the triangular lattice, and the prime on the sum restricts S to nearest neighbor vectors. We show 
an example for A = 0.25, h = 6.0 in figure 12. We have — | < 7^ < 1, which leads to the linear instability of the 
SF phase, which is denoted hn in equation(ll). Note that the mean field energies of the various phases lead to first 
order instabilities which make the linear instability irrelevant except at A = and A = i. 

We next turn to the transition line between the supersolid and Neel solid phases, denoted by h c x(A), a second order 
line.. Figure 13 shows the SW dispersion at the transition for A = 0.25, h — 0.9738. The soft density mode is now at 
k = and has a quadratic dispersion instead of the usual linear one. As we increase h and enter the Neel solid phase, 
the density mode becomes gapped, as figure 14 shows. The total energy (including spin wave corrections) in the Neel 
solid phase is independent of h, which demonstrates its incompressibility, and hence the exact filling of | throughout 
this phase. 



2. Kagome Lattice q = 

The Kagome lattice is a triangular Bravais lattice with a three element basis. If a is the nearest neighbor separation, 
then the Bravais lattice constant is a = 2a. The spin wave theory dispersions are quite similar to those of the 
triangular lattice, with the modes being softer (because of the lower coordination). Note that while the Bravais lattice 
is triangular, the symmetry is reduced (to C2,,) and the dispersion curves do not have zero slope at the zone edge 
(the X point). Figures 15, 16, and 17 show the dispersions in the superfluid, the supersolid, the Neel solid state, 
respectively. A noteworthy feature is the presence of a zero energy mode at h = along the TM direction in the 
supersolid phase. This mode disperses and has nonzero energy except along the TM direction. 



3. Kagome 9 Sublattice Structure 

The v^x V3 structure on the Kagome lattice is described by a triangular Bravais lattice of lattice constant 
a = 2\fi a with a nine element basis. The elementary lattice vectors are shown in figure 4. The spin- wave dispersions 
are unremarkable except for a flat mode whose energy vanishes as \fh in the small h limit. This is conected to the 
degeneracy mode, which is local up to harmonic order, and will be discussed in the next subsection. All the other 
features are quite similar to those of the triangular lattice as shown by figures 18, 19, and 20. 
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C. Gap of the Degeneracy Mode 



The degeneracy mode appears gapless in SWT (and in fact appears flat in the \/3 x -\/3 structure on the Kagome 
lattice), but acquires a gap to higher order in the SW expansion 25 38f| . We will first treat the triangular lattice, and 
go on to the case of the Kagome lattice. 



1. Triangular lattice 

Let us concentrate on the k = degeneracy mode of the triangular lattice. The easiest way to derive the gap is in 
the spin coherent state path integral language p2[ introduced in the previous section. 

Now imagine separating the path integral by first integrating all the modes except the k — degeneracy mode (call 
the angles 9 A and 4>°). 



Z = J D9 A D<f>°eiq>- J dr (iScos9° A ^- + \S 2 K{<p Q ) 2 ^ J J[D 



A'leM (3 2 ) 



fe^O 

where A' does not include the explicitly written terms involving <fio in the preceding factor. Since the classical ground 
state energy does not depend on 9 A , only the Berry phase term and a "potential" energy for cj> appear there (expanded 
to lowest leading order). If one neglects the integration over the remainder of the modes, the degeneracy mode appears 
to have a potential energy but no kinetic energy. This is analogous to a particle of infinite mass, and from the simple 
harmonic oscillator formula to = \J K/M, the oscillator energy is zero. This is why the mode appears gapless. 

However, as seen from figure 2, the integration over the rest of the modes (to leading order in the spin wave 
expansion) creates an effective potential which has a dependence on 6 A , rendering the mass M finite and the mode 
gapped. Since the effective potential is produced by SWT, it will be of order S, as opposed to the order S 2 potentials 
generated classically for the other modes. Thus the gap is order y/S, compared to the energies of order S seen in 
SWT. 

Now for the specific details. We concentrate on the neighborhood of 9 A , the optimal value of 9 A , and write 
9 A = 9* A + 69 A . It is easy to see that the 4>° that couples to this is <fi° — <fr A (k = 0) — (f> B (k = 0). The eigenvector 
corresponding to this mode is (</> A , 4> B , <p c ) = (§0°, — ^0°, 0). For the triangular lattice the classical "potential" energy 
of a (fP deformation yields 

A - = 3Wfl :( 1 + SSr) (33) 

Choosing the particular value A = 0.25 for illustration, we fit the curve in figure 2 near its minimum to obtain 
the term in the effective action i SC(59 A ) 2 , for which we obtain C ~ 6. Of course, a "potential" term will also be 
generated by the integration of the rest of the modes, but since it is order S, it can be neglected compared to the 
order S 2 term already present classically. We will choose the "coordinate" as Q — Ssin9 A (f> . The Berry phase term 
now looks like S s'm9' A 69 A (d(j> /dr), which enables us to identify P = 89 A as the conjugate momentum. We now write 
the action as 



S = I dr ( iP^- + \SCP 2 + \KQ 2 ) (34) 



,dQ 

dr 

where K = 3A(1 + | Sln 6 ° ) . We can now find the gap as the harmonic frequency of this oscillator: 



LO 



VSCK (35) 



The numerical value is uj ps 2.3-y/S. The Ka gome lattice structure introduces new considerations, which we now 
address. 
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2. Kagome Lattice 



The important difference that occurs for the Kagome lattice is that the effective action for the 9 A has a linear cusp 
at the minimum instead of a quadratic minimum, as shown in figures 21 and 22 for the q = and V3 x structures, 
respectively. This kind of minimum has previously been seen for the Kagome lattice in the fully antiferromagnetic 
Heisenberg case p| , ^2[ . We assume that the degeneracy mode is non-dispersing, which, at spin- wave order, is 
approximately true for the q = structure, and exactly true for the y/3 x V3 structure. This implies that the modes 
are local, and we assume this to be true even after quantum fluctuations have lifted them from zero energy. We can 
now interchange the roles of 9 and </> as P and Q and consider the quantum mechanics of the Hamiltonian 

^ + W\ (36) 

where A is of order 5. This leads to a gap in the degeneracy mode of order S% . A lifting of the degeneracy mode of 
order 5 a has previously been seen in in the Heisenberg case. 

Since the Kagome structure is always lower in energy, we will concentrate on it in the following, leaving 

a detailed analysis of the q = structure to a future publication. The degeneracy mode which was merely gapless at 
k = on the triangular lattice, but of nonzero energy at every other k, becomes completely flat at zero energy on the 
Kagome lattice with the V3 x V3 structure. The reason is fairly straightforward. Consider an ABABAB hexagon 
in figure 4b, isolated from the rest of the lattice by a ring of C sites. For the optimal angles it is easy to deduce 
that if 9 A = 9* A + 59, then 9 B — 9 A — 89 and that 9 C changes only to order S9 2 . Any coupling between the 59 A (r) of 
neighboring ABABAB hexagons must be mediated by the bordering 9 C , and must therefore be third order or higher 
in 59 A (r). This means that to quadratic order the 9 A fluctuations of the hexagon do not interact with the rest of 
the lattice. It is these local excitations that produce the flat mode. Note that even classically, the energy of this 
distortion is not zero if higher orders in 89 A are included. This is analogous to the flat mode in the q = structure 
of the isotropic Kagome antiferromagnct, where once again, the mode is flat only to harmonic order f35f| . 

Consider this degeneracy mode in the presence of a small field. It is easy to show that the field energy (per nine 
site unit cell) near the optimal angles is 

5£ fie i d = l h S 2 ^~_ A ^ cos 9 A (69) 2 . (37) 

We do not need to consider the bond energies of the spins to this order for the following reasons: 

• One can decompose the deviation from the optimal 9 A , 9 B , 9 C into S9 B — S9 BA + 59' B , and 59 c = 59 CA + 59' c , 
where the first part represents the change in 9 B and 9 C along the degeneracy direction, due to the change in 9 A , 
and the primed part corresponds to a change orthogonal to the degeneracy direction. 

• The optimal 59 B , 59' c in the presence of the field are of order h. This is because the bond energy is quadratic in 
59' B , 89' c while the field energy is linear in these quantities. Therefore, the energety difference due to the bonds 
will be of order h 2 , which can be neglected in comparison to the order h field energy at small fields. 

• There is no bond energy associated with a change along 59 A . 

The local 4> mode conjugate to this 9 A mode is (^5<p, — ^5<p, 0). The "potential" energy (per nine-site unit cell) of 
this mode can be calculated from the classical Hamiltonian exactly as in the triangular lattice case, and is 

U=\ S 2 A<fsin 2 9* A (2 + ^) (38) 

The Berry phase term for this unit cell, 35 sm9 A S9(d<j>° /dr), allows us to identify Q = 35 sin#*</>° as the coordinate, 
and P = 59 as its conjugate momentum. Then 



;(5, h, A) = oj(A)sVh 



This gap goes to zero as \fh in the limit h^O, and the mode becomes gapless in addition to being flat. 
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So far we have used the field to stiffen the degeneracy mode. Now consider the effect of quantum fluctuations at 
h = 0. We can find the effective potential to lowest leading order by computing the ground state energy for various 
6 A in SWT. This is plotted for A = 0.25 in figure 22. We assume that the mode will remain dispersionless, and 
therefore local, even when lifted from zero energy by quantum fluctuations. This assumption allows us to extract 
an effective Hamiltonian for each local degeneracy mode, which has been written down in cquation(35), where we 
identify the momentum as P = 3 S sin 9* A (jP and the conjugate coordinate as Q = 56, and where the numerical value 
of A is obtained from figure (for A = 0.25): 

AT 1 = 0.1944 

A = 0.025 S (40) 
We use a Gaussian trial wave function to obtain the approximate ground state energy of this Hamiltonian, yielding 

Of course, we expect this mode to disperse, but the calculation of the dispersion |^3|,[l4| is much harder than the one 
presented here, and will be pursued in future work. 



III. EFFECT OF QUANTUM FLUCTUATIONS ON ORDER PARAMETERS 



One of the motivations for this work has been to see if quantum fluctuations can disorder the system, creating a 
spin liquid, which would correspond to a ordinary liquid (with nonzero viscosity) for the bosons. To leading order in 
SWT, one can compute the average values of the spins as 

&) = «ata) - S)n MF . (42) 

The calculation sketched out in section II also produces the explicit Bogoliubov transformation, which can then be 
used to find the expectation values of bilinears. More explicitly, in terms of the matrix T which implements the 
transformation 

2K 

(e\{k)e j {k))= T L(k)T ja (k), (43) 

a=K+l 

where it is understood that i and j run from 1 to K (the number of sublattices), while a runs from K + 1 to 2K, the 
rank of T. Let us now turn to the different cases. 



A. Triangular Lattice 



Figure 3 shows a plot of the classical and quantum-corrected total S z and total S x , for h = 0. It is clear that the 
quantum-corrected value of the magnetization is very close to zero, independent of A. Exact diagonalizations of finite 
clusters for the anisotropic antiferromagnetic case A < Ej| suggest strongly that this is an exact statement. That 
is, the exact ground state at h = has zero longitudinal magnetization M z . An interesting fact about the triangular 
lattice allows us to map the A > problem (our model) to the A < model, which is the anisotropic antiferromagnet 
solved in SWT by KMFF. For very small A one can work to linear order in A, which means that one is working within 
the set of antiferromagnetic Ising model ground states. For nearest- neighbor spin flips, the Marshall sign property is 
obeyed by the wave- functions, arising from a partition of the set of ground states into disjoint even and odd states 
jh| . This means that for very small A there is evidence that the boson system would be exactly half-filled in the true 
ground state. However, to any higher order in A, or numerically for A not so small, the two models have no simple 
relationship with each other since the triangular lattice is not bipartite. 

Although total S x decreases as A decreases, it seems that superfluid order persists all the way, vanishing only when 
A = 0. In figure 23 we plot the quantum fluctuation corrections to the magnitudes of the spins on the three sublattices 
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(since 9 B = 6 C we plot only two values). As A— >0 the A sub-lattice remains stiff, while the B and C sublattice spins 
get reduced to about half their classical value (for S — j)- Thus, supersolid order on the triangular lattice survives 
quantum fluctuations at the spin wave level. 

Let us now turn to figure 24, which shows the classical and quantum order parameters as a function of A for h = 3. 
For A < i it is clear that though there are quantum fluctuation corrections to each of the spins, the total S z is not 
corrected, thus leaving the filling exactly |. For A > |, there are fluctuation corrections which do not affect the 
nature of the phase. 

Once again, we will concentrate on the lower energy \/3 x \/3 structure on the Kagome lattice, leaving an analysis 
of the q = structure to future work. 



B. Kagome 9 Sublattice Structure 



Figure 5 presents the naive results for the occupation numbers of the different sublattices in SWT for the \/3 x \/3 
structure at A = 0.25 as a function of h. Notice that as /i— >0 the fluctuation corrections diverge. This is a consequence 
of the flat mode, whose energy goes to zero as h—>0. In reality, as we have seen in the previous section, the energy 
fiat mode will be lifted in higher orders of SWT. While a quantitatively accurate analysis requires the computation of 
the dispersion of this degeneracy mode, we can make a rough estimate of the fluctuation corrections to the sublattice 
spins by assuming that the mode remains flat at the value u>i calculated in Section ILC.2. 

We need to determine the reduction in the magnitude of each sublattice spin due to quantum fluctuations. In terms 
of the equilibrium position 9* and the deviations from equilibrium 89 and 5(f> we can write the magnitude of the spin 
as 

p2 



\{S)\ = S{l-\(59f 

where we have once again used P = 3Ssm9 A (f. 
We know from the harmonic oscillator that 



\ sin 2 i 



2 ) 



1 



P 2 Q 2 \ 
2 18S 2 J 



(44) 



and Q — 9 A , and assumed a non-dispersing degeneracy mode. 



n = 

(Q 2 )o = 



2M- 
u 



2K 



and in the case h ^ 0, we find 



\(S)\=S- 



12SVh 



\ 



A(2 



sin 9^ 
sin 



cosf 



(2~A)2 
1-A 



h 
2 



cosf 



(2-A) 2 
1-A 



\ A(2 



sin 6* 



(45) 



(46) 



Clearly there is a divergence as h— >0, which is seen in figure 5. Of course, there will also be a contribution 
independent of h due to the other modes not associated with the degeneracy. 

Now consider the case h = 0. As shown in the previous section, the flat mode will be lifted from zero energy by 
quantum fluctuations and acquire an energy of order We will assume that the mode does not disperse (even 
though it will, weakly) in order to estimate the quantum fluctuations in the sublattice spins. We will also assume 
that none of the matrix elements change substantially. Thus the energy of the mode is the only significant factor. 
This means that we can estimate the contribution to the fluctuation correction to the spins by applying a field Hq 
such that the field-induced energy is the same as the energy induced by quantum fluctuations: 



A 2 



(47) 



Numerically, we find for A = 0.25 that /iq = 0.00057. This implies boson occupations of n A = n B = 1.66, and 
n c = 0.04. Since the boson occupations are larger than the value S = \, we find the fluctuations to be larger than 
the mean field value of the spin. Thus, despite the lifting of the flat mode due to quantum fluctuations, we find that 
fluctuations may be large enough to disorder the mean field ordering on the A and B sublattices. However, it is not 
certain that large fluctuations imply disorder. Furthermore, it is difficult to understand how the C-sublattice could 
remain stiff and ordered if the other two lattices are disordered. Strictly speaking, what this result tells us is that we 
have reached the limitations of spin-wave theory. 
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IV. CONCLUSIONS AND OPEN QUESTIONS 



We have found a broad region of parameter space where supersolid order is stable to leading order in mean field 
theory. The inclusion of spin wave corrections modifies this picture differently for the triangular and Kagome lattices. 
The supersolid is quite robust on the triangular lattice. Our interpretation of the supersolid corresponds with the 
ideas of Andreev and Lifshitz Q], with hopping vacancies undergoing Bose condensation. The role of lattice frustration 
is to prevent all the particles from condensing into a solid. 

We have also come across a general (though nonrigorous) argument which limits the search for the true ground 
state in a system with ground state degeneracy: Look for the points at which global conserved quantities commute 
in the Poisson Bracket sense with the generators of motion on the degeneracy subspace. 

On the Kagome lattice we find the \fi x \/3 structure to be more stable than the q = structure at all parameter 
values. Fluctuations seem much stronger here, and may even be able to destroy the long-range order assumed in mean 
field theory. This is a fertile region for numerical and experimental work. 

To what experimental systems might these considerations be applicable? One might be able to construct an array 
of Josephson junctions [jlTj which satisfy the conditions necessary for the existence of a supersolid. This means that 
the charging energy of each grain should be very high, and the nearest-neighbor charging energy should be higher 
than the Josephson coupling between the neighboring grains. Furthermore, only pair hopping should be relevant, 
which implies temperatures low compared to the bulk superconducting T c . Of course, since this is a two-dimensional 
system, the Bose condensate disappears for T^O, but power-law ODLRO is expected to remain. 

Another experimental system to which these results might be relevant is 4 He on graphite. A variety of orderings and 
transitions are known to occur as a function of temperature and coverage [Q. Also, steps in the superfluid density 
have been seen as a function of coverage (for multiple layers) for 4 He on graphite ]49| ] , and have been interpreted as 
resulting from correlation effects J50J . 

As it stands, this work is not applicable to the question of supersolidity in 4 He on a smooth substrate. In order to 
approach the continuum one would have to consider very low densities on the lattice, as well as long-range interactions 
(see ref |15 for an example of a continuum approach). 

We close with a number of important open questions. 

SWT seems sufficient for the triangular lattice, but not for the Kagome lattice. A formalism that can consider 
ordered and disordered states in a unified manner is necessary, perhaps a variant of the large- A approaches pl| , p6[ . 
Secondly, even for the triangular lattice, the question of vortices in the ground state is open. The spin-wave ground 
state has vortices, but they occur in nearest-neighbor pairs. A fruitful way to incorporate vortices might be to consider 
the effective theory of the gapless modes only, and include their interaction with vortices in a semiclassical manner. 

The original picture of ref |20|| has been confirmed by a strong-coupling expansion fl0|| . Strong-couplings expansions 
and variational approaches |52|j53|| are complimentary to the spin-wave expansion, and it would be interesting to 
investigate our model by these approaches. 

The effects of nonzero temperature on systems with ground state degeneracy can be very nontrivial. In particular, 
quantum and thermal ground state selection effects may compete to produce a sequence of phases and transitions as 
the temperature is raised |^9| . Clearly, thermal effects have to be elucidated before our results can be directly applied 
to an experimental situation. 

Finally, it is intriguing to harken back to the question raised by Anderson and Fazekas p4| and ask whether it is 
possible for a liquid state to be stable at zero temperature. 
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A. Appendix: Continuous Degeneracy not Arising from a Symmetry 



Consider a Hamiltonian function Ti.(x), where 
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X= (p 1 ,p 2 ,...,p N ,q 1 ,q 2 ,...,q N ) (48) 

is the vector of (dimensionless) coordinates and momenta. What does it mean that the ground state is continuously 
degenerate? Clearly, there must exist a one-parameter family of ground states, i.e. 

on 

dx 



= (49) 

f*/X(A) 



for all [i = 1, . . . ,27V, where A parametrizes the degeneracy submanifold, a curve in phase space defined by X(X). 
The degeneracy submanifold is one-dimensional in this example. Differentiating with respect to A gives 

i.e. the vector dX v /d\ is a null eigenvector of the Hessian matrix evaluated at the point X{\). 
As an example, consider the function 

H{r,<j>) = \(r 2 - l) 2 (a + focos0) (51) 

whose minimum is at r = 1 independent of <p. Note that </> is not cyclic in H, and its conjugate momentum (say r 2 , 
by definition) is not conserved as Noether's theorem does not apply. Nonetheless, there is a one-parameter family of 
ground states, parametrized by cj). The Hessian matrix, evaluated at r = 1, is 

( 2{a + bcos(f>) 0\ 

H -[ o oj 

and the vector ^ ( ) = ( 1 ) * S & nU ^ e ^S envector - 

Now consider the Gaussian fluctuations about the ground state. We are interested in the partition function Z = 
J V[p, q] e~ A , which is obtained from the action functional 

A = jTdr (ih 1> §7 + H & $J ( 53 ) 

where i = 1, . . . , N. We expand about one of the ground states, writing pi = Pi + Sp i} qi = Qi + Sqi, and, using 
summation convention, 

U = \8 Vl T tJ 5p 3 + | 5 Pl R.. 8q 3 + | 5q t R% S Pj + f 5q t V iS Sq 3 (54) 

with 

13 h\ d Pl d Pj ) p $ 13 h\ d Pl dq 3 Jpg v h\ dqi dq 3 J p g 

and R} is the transpose of R. Thus, the integrand of A is given by §fe M (H^ u + iJ^ u d T ) Sx„, where 

H=\ ] (56) 



and 



®NxN ^JVx AT 



J=\ \ ■ (57) 

V ^ NxN ®NxN 

The matrix H is brought to diagonal form by the basis change S: 
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(58) 

where 12 * = diag(wj t , . . . ,u) N ) is a diagonal matrix, and S^JS = J preserves the symplectic structure (commutation 
relations). Note that the diagonalized Hamiltonian is fl + p + q l q) rather than, say, \{p t p + q t Q 2 q). The two 
are, of course, related by a canonical transformation, provided the uj^ are nonzero.) The eigenvectors come in pairs 
| ip k ) which satisfy 

JH\^+) = -u>+\^) 



JH\^) =+^\4) 



as well as the orthonormalization condition 



_<^-|j|V,+ ) = (V^IJIVT) =*w 



U) 



k 



where k,l — 1, . . . , N. The N energy eigenvalues arc then given by e/t = hy uj^ 

Now consider the case where the degeneracy subspace is two-dimensional, parametrized by the coordinates (A, £). 
In our spin wave theory, we have that = fa and pi = 5(1 — cos #,), and the A invariance is realized as some 
relation among the colatitutes Oi(A) along the ground state submanifold, while the £ invariance is simply expressed 
as $i(£) = $i(0) + £. The matrix H (and hence JH) has two known null right eigenvectors - call them | N\ ) and 
| N 2 ) - corresponding, respectively, to the A and £ invariances. We have that 



dPi 



■ S sine, ^ 



iVO = I 1 = I (59) 



OQ, 



0, 



and 



OA / \ "N 



\ /n 
9« \ U N 



Note that the | N 2 ) eigenvector contains a 1 in each of its lower N entries. When we compute the overlap 
( N 2 | J | Ni ) , we obtain 

(N 2 \J\N l ) =-^9,^ = ^, (61) 



=i 



where M z is the total z component of the magnetization: M z = 5*cos0i. Unless the magnetization M z is 

extremized with respect to A, there is a finite overlap, and by appropriate normalization one can choose | Ni ) and 
| N 2 ) to be a | tjj ) pair. When the magnetization is extremized, however, the overlap is zero which means that 
these two null right eigenvectors belong to separate pairs, and there are at least two zero modes. 

Therefore, the two modes are independent only when dM z /d\ — 0, that is to say, when the magnetization M z is 
extremized. 

In the general case both the conserved quantity and the generator of motion along the degeneracy subspace may 
have components along both p- and g-subspaces. It is then easy to sec that the two vectors will have zero overlap only 
if the conserved quantity commutes (in the Poisson Bracket sense) with the generator of motion along the degeneracy 
subspace. 
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FIG. 1. Mean-held phase diagram for the triangular lattice. The Kagome lattice phase diagram differs only by a rescaling 
of h. Heavy lines denote first order transitions, light lines second-order transitions, and dashed lines denote linear instabilities. 
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FIG. 2. Ground state energy of the triangular lattice at A = 0.25 as a function of 9 A - The minimum is quadratic. 

FIG. 3. Classical and quantum magnetizations for the triangular lattice, S z and S x per site, as a function of A at the true 
ground state. The full quantum-corrected S z is always close to zero, while the quantum-corrected S x is always nonzero. 



FIG. 4. The two long range ordered structures on the Kagome lattice, (a) the q — structure and (b) the R3 structure. 

FIG. 5. Classical and quantum-corrected values of the magnetizations S z and S x per site as a function of h at A = 0.25 for 
the R3 structure on the Kagome lattice. Notice the divergence as h —> 0. This is an artifact of SWT. 

FIG. 6. Optimal 9 A = B and 6 C as a function of ft for A = 0.25. Note the continuous approach to collinearity at the 
supersolid-solid transition, and the discontinuous change of the angles at the solid-superfluid transition. 



FIG. 7. The Brillouin Zone of the triangular lattice and the reduced zone for the sublattice structure. 



FIG. 8. Triangular lattice SW dispersion for A = 0.75, h = 0. 

FIG. 9. Triangular lattice SW dispersion for A = 0.51, h = 0. Note the quadratic minimum which is nearly gapless. This is 
the linear instability that leads to the SS phase. 



FIG. 10. Triangular lattice SW dispersion in the SS phase at A = 0.25, h = 0. 



FIG. 11. Triangular lattice SW dispersion in the SS phase at A = 0.25, h — 0.5. The degeneracy mode has become gapped 



FIG. 12. Triangular lattice SW dispersion in the CS phase at A = 0.25, h = 6.0. Note the single gapless density mode. 

FIG. 13. Triangular lattice SW dispersion at the SS-MI transition at A = 0.25, h = 0.9738. Note the quadratic dispersion 
of the gapless density mode. 



FIG. 14. Triangular lattice SW dispersion in the MI phase at A = 0.25, h = 1.5. 



FIG. 15. Kagome q = SW dispersion in the CS phase at A = 0.25, h = 4. 

FIG. 16. Kagome q = SW dispersion in the SS phase at A = 0.25, h = 0. Note the flat chharacter of the mode along the 
TM direction. 



FIG. 17. Kagome q = SW dispersion in the MI phase at A = 0.25, h = 1. 



FIG. 18. Kagome R3 SW dispersion in the CS phase at A = 0.25, h = 4. 
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FIG. 19. Kagome R3 SW dispersion in the SS phase at A = 0.25, h = 0. Note the completely flat zero energy mode. 



FIG. 20. Kagome R3 SW dispersion in the MI phase at A = 0.25, h = 1. 



FIG. 21. Quantum energy correction versus 8 A for the Kagome q = structure. Note the cusp at the minimum. 

FIG. 22. Quantum energy correction versus A for the Kagome R3 structure. Once again the leading order potential has a 
linear cusp. 

FIG. 23. Quantum fluctuation corrections to the magnitudes of the spins on the A and C sublattices as a function of A at 
h — for the triangular lattice. 

FIG. 24. Classical and quantum values of S z and S x for the triangular lattice at ft = 3 as a function of A. Note that for 
A < i we are in the MI phase, and the total magnetization is uncorrected from its mean-field value of §. 
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